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ABSTRACT 


A significant area of interest in design of complex structures 
involves the study of multi disciplined problems. The coordination of 
several different intricate areas of study to obtain a particular design 
of a structure is a new and pressing area of research. 

In the past, each discipline would perform its task consecutively 
using the appropriate inputs from the other disciplines. This process 
usually required several time-consuming iterations to obtain a 
satisfactory design. The alternative pursued here is combining various 
participating disciplines and specified design requirements into a formal 
structural computer code. The main focus of this research is to develop 
a multidisciplines structural tailoring method for select composite 
structures and to demonstrate its application to specific areas. 

The development of an integrated computer program involves the 
coupling of three independent computer programs using an executive module. 
This module will be the foundation for integrating a structural optimizer, 

a composites analyzer and a thermal analyzer. 

With the completion of the executive module, the first step was 

taken toward the evolution of multidiscipline software in the field of 
composite mechanics. Through the use of an array of cases involving a 
variety of objective functions/constraints and thermal -mechanical load 
conditions, it became evident that simple composite structures can be 
designed to a combined loads environment. 
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CHAPTER I 


INTRODUCTION 


A major area of interest to date involves the study of 
multidisciplinary problems. Multi discipline is defined here as the 
coupling of technical disciplines such as structural analysis, composite 
mechanics, structural optimization, and heat transfer. It can be further 
described as the coordination of several different intricate areas of 
study to obtain a design of a particular structure which will concurrently 
satisfy these fields of interest. The study of fiber composite structures 
is the particular problem of concern here and will involve the four 
complex areas of study described above. The ability to have these 
disciplines interact with each other simultaneously on a specific 
structural design problem is a new and urgently needed area of research. 

In the past, each discipline would tailor its research to the 
particular interested field of study. If information from other related 
areas were needed, the task of collating the input from the relating field 
of study involved several time-consuming iterations to complete the task 
at hand. With the development of multidisciplinary structural tailoring 
methods into a formal structural computer code, a new tool for design of 
composite structures is within the grasp of the design engineer. 
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Though this research is limited to composite structures and a few 
of its related disciplines of study, it must be stated that there are 
potentially no limitations to the software development of integrated 
multidisciplinary research. 
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CHAPTER II 


RESEARCH OBJECTIVE AND APPROACH 

The main focus of this research is to develop a multidisciplinary 
structural tailoring method for select composite structures and to 
demonstrate its applications to specific cases. 

The underlying objective will be to develop an integrated computer 
program that will couple together three independent computer programs 
using an executive module. The software packages used have been developed 
at, or for, the Lewis Research Center in the recent past. The executive 
module developed here for the purposes of this research will be the 
foundation for integrating the three codes: STAEBL/GENCOM, ICAN and the 
thermal analyzer section of CSTEM. This integrated computer program will 
complement the conventional approach which presently requires substantial 
professional and computer time to obtain acceptable designs. 

The computer code will be subsequently used for the 
thermostructural tailoring of select composite structures. An array of 
cases which involves a variety of objective functions/constraints and 
thermal -mechanical loading conditions has been selected. 
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The array of cases includes four types of material systems, each 
with a fiber volume ratio of 0.60. The material systems are: 

1) AS/IMHS 

2) HMSF/IMHS 

3) HMSF/IMHS//SGLA/IMHS (intraply system) 

4) HMSF/IMHS and SGLA/IMHS (interply system) 

Definitions: 

AS: Graphite fiber 

HMSF: High modulus surface treated fiber 

SGLA: S-glass fiber 

IMHS: Intermediate modulus high strength matrix 

Intraply: 20% of each ply in layup will have the second system 

integrated within it 

Interply: 20% of plies in layup will be second system 

The optimizer within STAEBL/GENCOM uses design variables, decision 
variables and constraints to obtain optimal designs of these composite 
material systems. 

Design variables are perturbed at each iteration to give the 
optimizer a direction and magnitude of change toward a design which is 
considered to be optimal. Ply angles, ply thicknesses and composite 
thickness are design variables in these cases. 

Decision variables are used by the optimizer in STAEBL/GENCOM as 
the object to be optimized, this is commonly known as the objective 
function. The bending modulus of the composite, composite weight and 
composite displacements are typical decision variables. 
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The constraints are needed to define a feasible region in which 
the final design must satisfy. Modified distortion energy failure 
criteria of each material property will be used as the main constraints 
in the cases of concern. Another important constraint to be considered 
will be deflections of the structure. 

The type of thermal loads applied will be a temperature gradient 
along the length, across the width and through the thickness. The thermal 
analyzer is used to perform this task. Each gradient will be applied 
separately and in combination to each material system mentioned. 

The orchestration and coordination of input variables as well as 
variables passed between the codes will be explained further in the 
following chapters. Also, information on validation cases with 
predictable results will be collated and discussed. Finally, the 
demonstration cases involving the mentioned material systems will be run 
with the appropriate loading conditions. The results from these cases 
will be compared and analyzed. 
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CHAPTER III 


THEORY AND DESCRIPTION 


3.1 I CAN 

This section of the analyses will be used as the composites 
analyzer. This branch is necessary for analyzing the composite structure 
prescribed by the user in the initial design as well as all subsequent 
designs of the structure determined internally by the code. The composite 
structure is more precisely defined as a multilayered structure composed 
of unidirectional fiber- reinforced composite plies or layers. ICAN 
(Integrated Composites Analyzer) was predominantly designed to analyze the 
hygrothermomechanical response/properties of these multilayered 
composites. The types of layers recognized by the program are: 1) a 

standard composite system that consists entirely of a primary composite 
system made of one type of fiber and matrix, 2) an intraply hybrid 
composite system that is made up of a primary and a secondary composite 
system arranged in a prescribed manner “within" a ply, and finally, 3) an 
interply hybrid composite system consisting of different material systems 
(fiber and matrix) within the composite structure. 

By definition, a multilayered fiber reinforced composite structure 
is not homogeneous in the microscopic or macroscopic level, thus the 
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advantage to using I CAN comes to light. By employing the theories of 
composite micromechanics and laminate theory, this composites analyzer can 
portray the essentially heterogeneous composite structure as a synthetic 
homogeneous structure. The discussion of these theories involved in the 
analysis of a composite structure is expounded on in greater detail in the 
following sections. 

3.1.1 Composite Mechanics 

The equations used in ICAN that related the ply properties to 
the constituent (fiber and matrix) properties of the ply or lamina were 
derived through the use of composite micromechanics. This branch of 
composite mechanics was formally structured and based on certain 
assumptions pertaining to geometry of the structures and the fundamentals 
of solid mechanics. In the derivation of the equations, the main 
assumptions made were: 1) the ply resists in-plane loads as depicted 
schematically in figure 3-1, and 2) the ply and its constituents behave 
linearly elastic to the point of fracture as shown in figure 3-2. These 
equations for the material properties of a unidirectional fiber-reinforced 
lamina that were based on the constituent properties of the ply were 
proposed in simple form by Chamis (Ref. 1). 

Figure 3-1 shows the coordinate axis system used in the definition 
of the different material properties of the lamina. Note that the 
predicted equations will represent values for the equivalent homogeneous 
ply that is assumed to be transversely isotropic in the 2-3 plane. 
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Figure 3.1 : Typical Fiber Composite Geometry 
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Rgure 3 .2 : Stress -Strain Behavior of Unidirectional Fiber Composite 




Micromechanistic geometric relationships are described in figure 
3-3. Represented in these and all future equations: k denotes constituent 
volume fractions, p and X denote density and weight ratios respectively, 
the subscripts f, m, and l denote the affiliation with the fiber, matrix 
and ply respectively. The composite micromechanics equations for 
predicting mechanical, thermal and hygral properties are shown in figures 
3-4, 5 and 6 respectively. Also included in this study are the composite 
micromechanics equations for strength, and environmental effects of 
moisture and temperature on the material properties (Ref. 2). 

Another aspect of composite micromechanics important to this study 
is the stress failure criteria (Ref. 3). The importance is brought to 
light by the observation that a combined stress state will limit the 
strength of a ply much more than a simpler uniaxial load condition. One 
of the failure criteria output by ICAN is the modified distortion energy 
(MDEIE) shown in figure 3-7. The optimizer within STAEBL will determine 
the stress constraint that represents the maximum value for each ply that 
will most effect the design in an adverse way. It can be observed from 
figure 3-7 that the MDEIE stress failure criteria determines failure on 
a minimum value, not the maximum needed by the optimizer. Therefore, the 
following equation will be adapted from figure 3-7 and will be the 
representation value used within the optimizer for establishing the stress 
constraints of the composites. 

F « 1.0 - MDEIE 
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The optimizer determines the minimum value for "F" from each ply of the 
composite, then comparing that value to 1.0, the decision if a ply has 
failed is made. If "F" is less than 1.0 there is no failure, if "F" is 
equal to 1.0, the ply is on the verge of failing, and if "F“ is greater 
than 1.0, the ply has failed. 

3.1.2 Laminate Theory 

Laminate theory is the aspect of composite mechanics which will 
relate the individual lamina and its properties to the overall composite 
structure and its properties. In this discussion, the terms composite, 
laminate and composite structure will be used interchangeably, as will be 
the labels ply, lamina and layer. 

The stress-strain relations for a unidirectional ly reinforced 
lamina are contiguous to an orthotropic material under a plane stress 
state, or more simply stated, a 33 * a 23 * a 31 * 0. The derived stress- 
strain relations of the generalized orthotropic ply are defined with 
respect to the principal material axis of the ply, this relation then has 
to be transferred to the global structure coordinate axis of the 
composite. 

If each ply or lamina of the laminated composite structure is 
considered to be a thin plate, then according to Kirchoff's hypothesis for 
plates, the lamina is assumed to be in a plane stress state. Also, for 
the linear case of a deformed plate, the strains of each ply are related 
to the extension and curvature of the mid-surface (Ref. 4). The 
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Figure 3.3: Composite Micromechanics : Geometric Relationships 
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Rgure 3.4: Composite Microtnechanics : Mechanical Properties 
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Figure 3.5: Composite Micromechanics : Thermal Properties 
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Figure 3.6: Composite Mteromechanics : Hygral Properties 
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Figure 3. 7 : Modified Distortion Energy (MDEIE) Failure Criteria 


relationship between the individual ply strains and that of the composite 
is shown in figure 3.8 where z k is the perpendicular distance from the 
laminate mid-plane to the mid-surface of the k th ply. Even though the 
strain variation is linear, the stress variation through the thickness of 
the composite might not be since the stress-strain relation can be 
different for each ply. 

It can be seen from figure 3.8 that the resultant stresses and 
moments acting on the composite can be derived. This is done by 
integrating the individual ply stresses over the thickness of the laminate 
(Ref. 4). 

This brief review of laminate theory, or better stated, classical 
laminated plate theory, should be understood to be a superficial study of 
this intricate field. For a more integral study of the important elements 
that comprise laminate theory, the ICAN User's Manual (Ref. 5) cites many 
references. 

In summary, figure 3-9 will state the important equations for 
laminate theory. These equations will include the hygral and thermal 
relations as well as the mechanical relations. 

3.2 Thermal Analyzer 

The thermal analyzer used here was developed using the 
fundamentals of heat transfer (Ref. 6). Heat transfer falls in the 
general classification of a transport phenomenon. Also included in this 
broad classification are the areas of mass transfer, momentum transfer or 
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Figure 3.9: Summary of Governing Equations of Laminate Theory 



fluid friction and electrical conduction. The rate equation of all these 
unique transfer systems is very similar in that the flux is proportional 
to a potential difference. The flux is defined as the quantity 
transferred per area per time; and the potential difference, in the case 
of heat transfer by conduction and convection, is the temperature 
difference. From this, the definition of heat transfer is described as 
the transmission of energy from one region to another primarily as a 
result of temperature difference. 

Generally, heat transfer is split into three different modes 
which are referred to as conduction, convection and radiation. Of the 
three modes mentioned, convection is the only one that does not comply 
with the definition described above since mass transport is also 
involved, but "heat transfer by convection” is a widely accepted 
term. 

The thermal analyzer code has the ability to analyze each of the 
mentioned modes of heat transfer: conduction, heat transfer by convection 
and radiation. This code was adopted for use in this research due to its 
flexibility. Some of the features incorporated within the code involve 
types of heat transfer analysis (linear steady-state, nonlinear steady- 
state, linear transient and nonlinear transient), material types 
(isotropic and orthotropic), element types (8, 16 and 20 nodes) and units 
(metric or English). Additional details about what components were used 
and why will be explained in much greater detail in the next chapter. 
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3.2.1 Conduction 

Conduction is the process by which heat is transferred from a 
region of higher temperature to a region of lower temperature within a 
medium. That particular medium can be described as a solid, liquid or 
gas. Also, heat transferred between different media in direct physical 
contact where there is no heat transfer by mass movement can be considered 
conduction. The energy is transmitted by direct intermolecular collision 
without any appreciable displacement of the molecules. 

The internal energy of a system is the energy possessed by an 
element of matter due to the velocity and relative position of the 
molecules. This energy is directly proportional to temperature of the 
element according to kinetic theory. From this it is deduced that the 
greater the temperature and thus, the internal energy, the more rapidly 
the molecules are vibrating. If there is a difference in temperature 
between two adjacent regions, the molecules of greater energy will lose 
part of their energy to the molecules of lower energy. In fluids, this 
energy loss is by elastic collision while in metals the loss comes in the 
form of diffusion of electrons. 

Thermal conduction can be expressed by the following relationship: 

q 

— = -k V T 
A 

q = heat flux 
k * thermal conductivity 
V = vector gradient operation 
A = area 
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This theory was proposed by J. B. Fourier in 1822 and is a 
generalization of the experimental work done by 0. B. Biot in 1804 and 
1816. 


3.2.2 Convection 

Heat transfer by convection is the method of energy transport 
by the combined action of conduction, energy storage and mixing motion. 
Heat transfer between a solid surface and a fluid can best be described 
using convection, since this method of energy transfer involves mixing and 
diffusion as well as conduction. To facilitate in understanding this 
intricate mode of energy transfer, examine the situation of heat transfer 
to a fluid flowing inside a pipe. Three different flow regions are 
present for a fast flowing fluid. The laminar sublayer adjacent to the 
pipe wall uses heat transfer by thermal conduction. The transition region 
outside the laminar sublayer has eddy mixing as well as thermal 
conduction. Finally, toward the center of the pipe, eddy mixing is the 
prevailing method of energy transfer. 

When modeling heat transfer by convection, several steps must be 
considered. Assuming a surface has a temperature greater than the 
surrounding fluid, the first step is the energy transfer to the fluid 
particles in direct contact with the surface. This results in increased 
temperature of these fluid particles. These more active fluid particles 
then will move toward and mix with the lower temperature fluid particles. 
This fluid movement has a two fold effect; more lower temperature 
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particles can now be in contact with the surface of higher temperature, 
and the more energetic fluid will transfer some of its energy to other 
fluid particles. This transfer is a result of mass motion. At this time 
it can be pointed out that due to this resulting mass motion, this mode 
of transfer does not strictly depend on a temperature difference and, as 
stated before, does not conform to the strict definition of heat 
transfer. The overall heat transfer effect is in the direction of the 
temperature gradient and therefore is loosely classified as a type of heat 
transfer. The class of heat transfer by convection is dependent upon the 
manner of motivation. The one described above is considered free 
convection. The mixing motion is due to density differences caused by a 
temperature gradient. Another cause of mixing motion is due to an 
external effect. An example of this could be from the use of a pump or 
blower; this form of convection is called forced convection. 

In 1701, Sir Isaac Newton successfully analyzed this intricate 
and unique field of convection. He proposed the general Newton Rate 
Equation: 

q ' hA <T„ - T„ ) 

q = heat flux 
h = convection coefficient 
T w = surface temperature 
T. = fluid temperature 
A = area 
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Convection coefficient, h, was a consideration of fluid motion, fluid 
conductivity and the role of turbulent eddies. 

3.2.3 Radiation 

The final method of heat transfer is by radiation. This occurs 
when a high temperature body transfers energy to a low temperature body 
when the two are separated in space. Radiant heat transfer acts from the 
internal energy at a source first being changed to electromagnetic 
energy. This energy is then transmitted to a receiving material and 
changed into internal energy in that material. Unlike conduction and 
convection, radiation is not dependent upon a material to act as a medium 
to convey the energy between the two regions, actually thermal radiation 
is impeded by the presence of any intermediate material. Radiant heat 
emitted by a body is in the form of finite groups of energy called 
quanta. These quanta have motions in space similar to light propagation 
and therefore can be described in wave theory. 

The expression that describes thermal radiation heat transfer has 
the form: 

— = a T 4 
A 

q * heat flux 
T = temperature (absolute) 
a - constant 
A = area 
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In 1879, Stefan obtained this expression from empirical data of 
Timdall, and then Boltzmann, in 1884, used classic thermodynamics to 
derive it. The expression is more commonly known as the Stefan-Boltzmann 
Law. 


3.2.4 Finite Elem ent Formulation 

As stated previously, this module can analyze linear steady- 
state, nonlinear steady-state, linear transient and nonlinear transient 
conditions. The details disclosed next should give insight to the finite 
element formulation used in the thermal analyzer code (Ref. 9). The use 
of three-dimensional isoparametric solid elements will allow the 
temperature, T, within the element to be defined at time, t, in terms of 
nodal temperatures. 


T(x, y, z, t) * I H,(x, y, z) T,(t) 

where H, is the shape function and T, is the temperature at node, i. 

The general three-dimensional heat transfer condition in a body 
assumes the material obeys Fourier's law of heat conduction. The 
expression for heat flow equilibrium in the interior of a body can 
follow: 


8 8T 


—(K — ) 

8x ox 


+ 



8T 

ly 


+ 


8 3T 


—(K — ) 
8z 3z 


-9b 


(3-1) 
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where q b is the rate of heat generated per unit volume, the thermal 
conductivities of the principal axes are denoted by k x , k y , k z and T is the 
temperature of the body. On the surface of the body, the following 
conditions must be satisfied: 


T - T. (3-2) 

SI 


<JT 



where T e is the environmental temperature and q $ is the heat flow input to 
the surface of the body. The details of these boundary conditions are 
expounded upon in the user's manual for the thermal analyzer (Ref. 9). 

To solve the differential equation, 3-1, with the given boundary 
conditions, 3-2 and 3-3, a finite element scheme was developed. A 
variational formulation of the heat transfer problem was engaged. This 
formulation defined a function such that when stationarity of that 
function is invoked, the governing differential equations 3-1, 3-2 and 3-3 
were obtained. 

The user's manual (Ref. 9) outlines a general solution scheme for 
linear and non-linear, steady-state and transient problems. The 
development of incremental equilibrium equations for each problem was 
done using the modified Newton-Raphson iteration for heat flow 
equilibrium. 
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At each iteration, the stiffness matrix will remain constant while any 
change will be realized by the load vectors. 

The cursory description of the finite element formulation outlined 
above should be deduced to be a review/summary. The details are 
delineated in the thermal analyzer user's manual (Ref. 9) and can be 
referenced if an in-depth study is required. 

3.3 STAEBL : 

Structural tailoring of engine blades/general composite structures 
(STAEBL/GENCOM) is an extension and modification of the program STAEBL 
(Ref. 7). Most of the modifications were carried out in order to apply 
the program to general composite structures. Since the overall program 
logic follows that of the original STAEBL program, this part of the module 
will be described first followed by the capabilities of the GENCOM 
version. 

The STAEBL computer program was developed to perform engine fan 
and compressor blade optimizations through the use of realistic blade 
design constraints such as blade stresses, vibratory response, flutter and 
foreign object damage. The optimization is reached by tuning one to 
twenty design variables. Included are airfoil thickness at several 
locations, blade chord and certain blade internal configuration 
variables. 

Three component parts are required to perform a blade 
optimization: an optimization algorithm, approximate analysis procedures 
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for objective function and constraint evaluation and refined analysis 
procedures for design validation. 

Control Program for Engineering Synthesis/Constrained 
Minimization (COPES/CONMIN) optimization package is used as the 
optimization algorithm in STAEBL. This is a proven tool for optimizations 
with a small to medium number of design variables. The optimization 
method using COPES/CONMIN will be discussed in greater detail in the next 
section. 

The approximate analysis utilized is an efficient, coarse mesh, 
plate finite element procedure in STAEBL. How this approximate analysis 
was decided upon will be reviewed in greater detail in a later section of 
this chapter. This analysis can also provide blade frequencies, mode 
shapes, stresses under loads, flutter and foreign object damage, which are 
utilized to evaluate respective design requirements in the form of 
constraints. 

A refined analysis should be applied to the optimal design to 
assure that all constraints are satisfied. This task is not performed 
by STAEBL automatically, but is the user's responsibility by incorporating 
an existing design/analysis system. If the first choice of optimal design 
is found to violate one or more of the constraint values, the allowable 
for the constraint values must be adjusted to take into consideration the 
difference between the approximate and refined analysis. During the 
development of STAEBL, a fully satisfactory design was found for all cases 
studied. 
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To use STAEBL, the coordinate description of an initial design 
must be known. From here the optimization system will adjust the design 
variables to try to establish an optimal design. In addition to geometry 
variables, there are additional construction variables available for 
composite blades, some of which are material thickness and ply angle 
orientation. 

3.3.1 STAEBL/GENCOM 

The GENCOM version of STAEBL (Ref. 7) had several additions and 
modifications in order to incorporate composite structural analysis and 
also allow the user certain options not in STAEBL. The user now can input 
a NASTRAN blade geometry description and a second option replaces calls 
to the external math library, IMSL, with calls to subroutines added to the 
source code. 

Since composite blades can experience high temperatures and large 
temperature gradients, through-the-thickness as well as in-plane, a 
thermal stress analysis capability was added. The capacity to analyze 
thermal gradients will allow the composite properties to be temperature 
dependent. A similar adjunct involving moisture gradients was included 
in the new version of STAEBL. 

The composite construction variables added to STAEBL involve 
simple composite blade preprocessing. Up to seven material systems can 
be analyzed with the thicknesses and fiber angles of each material system 
incorporated as design variables. The design thickness is built up from 
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the outside toward the core with the thickness of each material system. 
Thickness of the core material system is adjusted to accommodate the 
balance of the design thickness. 

3.3.2 COPES/CONMIN 

As stated earlier, the COPES/CONMIN (Ref. 7) optimization 
package is used as the optimization algorithm of STAEBL. It was 
selected because of its versatility for structural optimization problems. 
It is applicable to both constrained and unconstrained minimization 
problems. 

The engineering design problem involved in STAEBL/6ENC0M is the 
determination of design variables which minimize a design quantity while 
satisfying a set of auxiliary conditions. This is generally a 
constrained minimization problem and is mathematically expressed as 
follows: 


minimize f(x) 


(3-4) 


where f(&) is the scalar function to be minimized. This function, f(&), 
is the objective function subjected to m number of auxiliary conditions 
known as the inequality constraints: 

9i(i) < 0, i « 1 m (3-5) 
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The vector of n number of design variables is represented by the 

quantity * = (x x n ) . Each design variable, x,, has imposed upon them, 

upper and lower bounds referred to as the side constraints. These side 
constraints can be represented in the form: 

L, < x, < U,, i * l,...,n (3-6) 

A feasible region is the combination of all feasible points in a 
design space. Each feasible point is any choice of design variables, x, 
which satisfy all the constraints, equations 3-5 and 6. 

The feasible region is bounded by a constraint surface or a locus 
of points which satisfy the equation g, (*) * 0 for some i. On one side 
of the surface, g^x) < 0, the constraint is satisfied, on the other side, 
9i(x) > 0, the constraint is violated. Any feasible points on the 
boundary, g, (*) = 0 are called bound points and that particular constraint 
is said to be active. Any points inside the feasible region are known as 
free points. If the objective function is at a minimum and the design 
point is in the feasible region, the solution to the problem presented in 
equations 3-4, 5 and 6 is known as an optimal feasible point. A feasible 
region is disjoint if it is made up of two or more distinct sets of 
feasible points. This could occur in any nonlinear minimization problem, 
that is, when f (x) or any g,(x)'s are nonlinear. In this instance, there 
can be multiple local minima and the global minimum is now the optimal 
feasible point. A structural optimization problem almost always finds a 
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solution on the boundary of the feasible region. In fact, it is usually 
at the intersection of two or more constraint surfaces and therefore the 
solution will have two or more active constraints. 

STAEBL uses the direct method, versus indirect method, to 
solve the constrained optimization program. In this method, the 
objective function and the constraints are evaluated independently and 
the constraints are treated as a limiting surface. Zoutendijk's method 
of feasible directions is an example of a direct method of solution 
and is incorporated in the exact analysis available for a solution 
technique. 

3.3.3 Finite Element Analysis 

STAEBL (Ref. 7) uses an approximate analysis procedure to 
establish an optimal design. This type of analysis is required to obtain, 
with reasonable accuracy, a design as quickly and as efficiently as 
possible. At this point, the acquired design is then evaluated using a 
more rigorous refined analysis for comparison with the approximate 
analysis. This comparison will allow the user to become aware if the 
optimal design is valid. 

NASTRAN finite element analysis is used for the refined analysis, 
and a beam analysis procedure was initially used for the approximate 
analysis. It was found that to properly model complex blades, a plate 
finite element was required for the approximate analysis. It can be shown 
by using the same plate finite element analysis as NASTRAN, with a 
comparatively coarse mesh, satisfactory results were obtained with 
competitive run times to the original beam analysis procedure. 
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Therefore, for the approximate analysis procedure and using 
NASTRAN plate technology, a self-contained plate finite element analysis 
procedure was developed. NASTRAN technology was chosen for several 
reasons: 1) proven computational efficiency, 2) established successful 
correlations with test cases, 3) the convenience of the input and output, 
and importantly, 4) the compatibility with NASTRAN refined analysis 
procedures. 

With the usage of a plate bending triangle, the approximate 
analysis in STAEBL preserves the similarity with NASTRAN. This 
isoparametric element of the QUAD4 family is very similar to the NASTRAN 
TRIA3 element. This TRIA3 element is a reduced integration triangular 
plate bending element with the following features: 1) recognition of 
thickness taper, 2) properly stacked triangular plate element meshes to 
simulate airfoil pretwist and camber, 3) composite material capabilities 
(using lamination theory), 4) element differential stiffness, and 5) 
lumped masses are employed, assuring a diagonal mass matrix, for storage 
efficiency. 

Since the approximate plate finite element analysis procedure is 
derived from NASTRAN, it uses NASTRAN-format bulk data for its input and 
produces NASTRAN-format displacement and stress output. Duplicate pre- 
and postprocessing can be used for both the approximate and refined 
analysis; however, since all storage and processing of the approximate 
analysis take place in the core, the program will have limited capacity, 
but the allowable size can be shown to be sufficient for relative 
accuracy. 
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The finite element mesh and final Guyan reduction pattern is shown 
on a model in figure 3-10. With this comparatively coarse mesh in 
combination with the Guyan reduction, it reduces a 330 degree-of-freedom 
model down to 24 degrees-of -freedom. This Guyan reduction procedure 
(Ref. 8) has been shown to reduce the number of degrees-of-freedom with 
minimal loss of accuracy in a dynamic analysis. In STAEBL, the reduced 
or omitted degrees-of-freedom, u 0 , and the remaining or analysis degrees- 
of-freedom, u 4 , are related to static loads according to: 



K.o 


“a 







* 


K oa 

K oo 


u 0 


F o 


neglecting the forces F e , 

K ] « [GJ [uj 


where 

[GJ - - [KJ - 1 [KJ 

the reduced stiffness matrix thus becomes: 

[*J = [*J + [*J KoJ 

The reduced mass matrix is determined by equating the kinetic 
energies before and after the reduction and becomes: 

[M„] - [MJ ♦ [MJ [GJ + [GJ 1 (M„ ♦ GJ 
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Figure 3. 10 : STAEBL Guyan Reduction Pattern & Boundary Conditions 











After the stiffness and mass matrices have been reduced, they are, 
in general, symmetric but full, yet small in size, and the unsymmetric 
eigenvalue problem becomes: 

-« 2 { u, ) + [HJ- 1 [KJ j u, ) - (0) 

Both eigenvalues and eigenvectors are found for the reduced size problem. 
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CHAPTER IV 

COMPUTER IMPLEMENTATION 
4.1 General Framework 

The essence of this section will be the discussion of each module 
of the program, how it was used, where it was used and why it was used. 
Also addressed herein will be any modifications or adaptations to each of 
the modules. As discussed earlier, the program consists of four parts: 
the executive module, ICAN, thermal analyzer and STAEBL/GENCOM. The 
theories and backgrounds of the last three modules were discussed in the 
previous chapter. This segment will examine closely the interaction that 
each module has with each other as well as its collaboration with the 
executive module. 

Figure 4-1 displays the general framework and flow of the 
program. To start, all variables for an initial design are input, this 
initial design sets up the geometry and design criteria for ICAN, 
STAEBL/GENCOM and the thermal analyzer. Included are connectivities, 
boundary conditions and mechanical loads applied to the structure. As 
earlier specified, STAEBL/GENCOM will use NASTRAN-type formats and figure 
4-2 shows the arrangement of nodes and plate elements used in this study. 
The thermal analyzer will use an 8-node brick element, figure 4-3 displays 
the geometry of the structure as viewed by this module. 
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Figure 4. 1 : Flowchart 


38 












39 


Figure 4. 2 : STAEBL Finite Element Mesh 
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Figure 4. 3 : Thermal Analyzer Finite Element Mesh 





After the input variables are read, the next steps of the program 
involve ICAN and the thermal analyzer. ICAN is called to provide the 
thermal conductivities of the composite material at room temperature. 
These thermal conductivities are employed by the thermal analyzer to 
generate the temperature profile in the structure. This temperature 
profile is dependent upon the thermal loads applied to the structure. 
ICAN is now engaged to produce the element membrane and bending properties 
along with the failure criteria and thermal stresses. The location of 
this latest call to ICAN is dependent on the type of thermal load 
condition encountered. The call will be denoted with dashed lines on 
Figure 4.1 between ICAN and the other two modules. If a thermal gradient 
exists along the length or width of the structure, each plate element may 
have a different temperature gradient. In view of this fact, combined 
with the reality that many of the properties are temperature dependent, 
they must therefore be calculated for each element and so this final call 
to ICAN comes from STAEBL/6ENC0M. If it is found that the temperature 
gradient through the thickness is the same for each element, then one call 
from the thermal analyzer is all that is needed to generate the properties 
for the entire structure. 

The subsequent stage of the program involves the optimizer and 
finite element analysis of STAEBL/GENCOM. The optimizer will perturb each 
design variable from the assigned initial design, then the analysis 
outlined in the previous paragraph is applied to each of these designs. 
The finite element analysis uses the calculated properties to output the 
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displacements for each node and the stresses for each element. These 
elemental stresses are converted to internal mechanical loads to be used 
in the final ICAN analysis for determining failure criteria. By comparing 
the output of each design, a gradient is established that will move the 
global design toward an optimum as determined by the assigned decision 
variable and constraints. 

Once a global design is decided on, it is then compared to the 
previous global design. If there is no relative change in the design 
variables, the global design is said to be optimal and all pertinent 
information is output. If the design is determined not to be optimal, the 
process starts over and continues until an optimal global design is 
established. 

4.2 ICAN 

The only change in ICAN apparent to the user is the load cards. 
Originally there were three "PLOAD" cards for each loading condition. The 
first card contained entries for the membrane loads along with their 
orientation. The second card contained the bending resultants, and the 
last card, the transverse shear resultants and the transverse pressures. 
These cards no longer are needed considering the loads now are being 
generated from the stresses output in STAEBL/GENCOM. In place of the 
three "PLOAD" cards a single "ELOAD” card will be used and will have the 
format shown below. 


1 8 

9 16 

17 24 

25 32 

33 40 

MNEMONIC 

NSTR 

MEND 

NIHC 

NANL 

ELOAD 

9 

Ifi 

Z 

1 
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The card in this group starts with the mnemonic "ELOAD". Each 
entry in the card will be defined as follows: 

NSTR: starting element number 

NEND: ending element number 

NINC: increment of element numbers between the NSTR and NEND 

NANL: type of thermal gradient expected: 

1: gradient through the thickness only 

2: gradient along the length only 

3: gradient along the width only 

4: any gradient combination of 1, 2 or 3 

Note that the element numbering system used here is the same used in 
STAEBL/GENCOM, see figure 4-1. 

Outwardly, the use of the "ELOAD" card is a minor change, but in 
reality it has multiple effects. First, it informs ICAN which elemental 
mechanical loads to apply to its analyses since the stresses within each 
element may differ depending on the type of global load conditions 
applied. Secondly, the thermal analyzer uses the gradient information 
when assigning temperature output to each layer in each element. Finally, 
STAEBL/GENCOM utilizes the data in determining the pertinent failure 
criteria to be conveyed to the optimizer for possible constraint 
evaluations. 

4.3 Thermal Analyzer 

The sole change discernible in the thermal analyzer was the method 
by which the thermal properties are input. Initially, these properties 
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were contained in the thermal analyzer input deck. They are now being 
generated in the ICAN module and passed to the thermal analyzer. This 
procedure is considered valid only because the composite analysis done by 
ICAN will render the essentially heterogeneous composite structure intp 
a quasi -homogeneous structure where the thermal properties along the 
structural axes are known. The thermal analyzer uses these properties to 
evaluate the temperature profile within the structure due to the thermal 
loads applied. This temperature profile will be utilized by the other 
modules in the structural analysis procedure. 

4.4 STAEBL/GENCOM 

The modifications to STAEBL/GENCOM were much more extensive than 
the other two modules. All but two of the input files were either 
eliminated or hard coded into the module. At present, the needed input 
data includes any optimization information, geometry and connectivities 
of the initial design, along with the initial layup and design 
information. 

The boundary conditions are needed as part of the STAEBL/GENCOM 
input file. For purposes of this research, nodes 1-5 (see figure 3-10) 
are restricted from movement in all directions except along the width. 
This will allow for thermal expansion of these nodes. The boundary 
conditions that involve the Guyan reduction pattern (see figure 3-10) are 
entered as ASET card information directly into the code. 

The mechanical loads applied to the structure are input through 
the SAVX file. The loads are entered using English notation and will be 
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applied in much the same manner as the boundary conditions. There will 
be 55 lines of input to this file, one for each node. Each line has six 
entries. The first three correspond to an extension load applied in the 
directions of the main axis, that being along the length, through the 
thickness, and along the width. Figure 4-2 shows the positive directions 
of each. The final three entries correspond to moment loads applied in 
the same direction as the extension loads. These mechanical loads 
combined with thermal and hygral loads are employed by the finite element 
within STAEBL/GENCOM to output the stresses (a xx , a yy , o xy ) at the top and 
bottom of each element. Owing to the fact that this is a linear case 
study, a combined stress formula can be incorporated with this stress 
output to calculate the loads needed for the composites failure analysis. 

The combined stress formula uses the coalition of stresses due to 
axial and bending loads. Therefore, 


a 


N 

A 


± 


M c 
I 


N = 
M = 
A * 
I = 

c = 


top or bottom stress, i.e. o T and o B 

axial loads (N c „, N cyy , N cxy ) 

bending loads (M cxx , M cyy , M cxy ) 

area per unit length, i.e., thickness (t) 


moment of inertia per unit length. 


t 3 

i.e., — 
12 


t 

distance to stresses, i.e. — for top and bottom 

2 
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and so: 


6 M c 


(4-1) 


N c 


6 M, 


(4-2) 


by adding equations 4-1 and 4-2: 


°c + °c 


B 2 


It follows that: 


- K T ♦ o c B ) 


(4-3) 


by subtracting equations 4-1 and 4-2: 


■ 


12 M c 


it then can be shown that: 


t 2 t B 

M c * (o c T + o e B ) 

12 


(4-4) 


Equations 4-3 and 4-4 are used for the calculation of N cxx , N cyy , 
N exy and M cxx , M cyy , M cxy at each element. These membrane loads are 
transferred to the ICAN module for the composite analysis of any element 
in question. The elements in consideration have been defined on the 
"ELOAD" card located in the ICAN input deck. 
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Finally, in addition to the above changes, alterations were made 
to the global variable array. This array contains all the possible 
decision variables, design variables and constraints used by the optimizer 
in determining an optimal design. The new variables include the NASTRAN 
bending equivalent elastic coefficients (Gil, G22, G12), maximum static 
nodal displacement in the directions of the main axes (x, y, z), overall 
composite thickness and the combined stress failure criteria of each 
composite material. 
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CHAPTER V 


CASE STUDY APPLICATION 

In the previous chapter, the general framework of the executive 
module was described. Also included was an explicit explanation of how 
each of the other modules were incorporated into the program. The next 
step will be the verification of the integrated program by applying a set 
of test cases. This chapter will characterize the design factors involved 
in this validation process. It will also define any involvement these 
design factors may have with the predescribed demonstration cases. It 
would be appropriate at this time to state that the test cases described 
here will only be a verification that the executive module has 
successfully integrated all the other modules involved in this research. 
It is assumed that the modules have been individually verified. 

The typical cross-section of the composite structure involved with 
each test case, and later the demonstration cases, is shown with figure 
5.1. This figure shows the properties of the cross-section that were used 
as design variables: composite thickness, material thickness and ply 
angle. This generic section may be adjusted considering the allowance of 
these design variables to change over a predetermined range. The range 
is dependent upon which property of the cross-section is being addressed. 
A maximum composite thickness of 0.625 inches was used for all cases 
studied. For the individual material systems, the maximum thickness was 
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Figure 5.1 : Typical Composite Section 
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COMPOSITE THICKNESS 






0.050 inches. Material system No. 7 was the only exception, a thickness 
of 0.025 inches was applied. The material systems are built up from 
individual plies, each 0.005 inches thick; the ply angles can vary from 
+90.0 degrees to -90.0 degrees. If a particular material system is not 
needed due to the relationship of the overall composite thickness to the 
combined thicknesses of the seven material systems, the material system 
number in question will be labeled "not used" in the tables referred to 
later in this chapter and the Appendix. This is done in order to present 
the results in a consistent format. 

Several load conditions will be attached to this structure, and 
for each load condition an appropriate objective function is chosen. The 
objective function for the test cases, as well as the demonstration cases, 
will be a bending modulus that is considered to best model the equivalent 
of a NASTRAN-obtained bending modulus for the entire structure. The 
optimal value of the bending modulus in question will be the largest 
possible value. 

The choice of design variable is dependent on the objective 
function used; a design variable that has a noticeable effect on the 
objective function is most desirable. The ply angle direction of each 
material system was the first choice of a design variable because of its 
strong effect on the bending modulus of a composite structure as well as 
the failure criteria of the individual plies. The ply angles of each 
material system will be measured off the 1-1 (x-x) axis for the validation 
cases as well as each demonstration case. When an attempt was made to 
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facilitate a clear design, it was found that ply angle direction was not 
as effective as first thought in some of the demonstration cases. 
Material system thickness and the overall thickness of the composite were 
some of the acceptable substitutes used. For the purpose of program 
validation, the ply angle directions were determined to be suitable design 
variables. The use and results of the other design variables mentioned 
will be expounded upon in the section dedicated to the demonstration 
cases. 

The final factor involved in controlling the design tailoring of 
a composite structure will be choosing the constraint variables. For the 
test cases, the combined stress failure criteria and maximum static 
deflection of the structure will be used. With regard to the 
demonstration cases, only the combined stress failure criteria was deemed 
relevant for design purposes. The stress failure criteria is based on the 
failure of an individual ply due to an applied stress. This translates 
into a failure of the entire structure. 

5.1 Validation Cases 

Tables 5.1 through 5.5 (pages 52 through 56) depict the form of 
each validation (test) case to be presented. This same descriptive format 
is used to document the demonstration cases. The top of each table shows 
the fiber/matrix system as well as the load conditions applied to the 
composite structure. The balance of each table will display the initial 
and final design values for the objective function, the design variables 
and the constraint variables. 
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ONSTRAINT | OESjGN I OBJECTIVE 


TABLE 5. 1 : TEST CASE # 1 


FIBER /MATRIX SYSTEM 

LOAD CONDITION q a 

*< 

r 

X 


T300 / IMHS 
FVR: 0.62 

TEMP. EACH NODE £ * 
T U =250.0°F M 

s 

OPTIMIZATION VARIABLES 

INITIAL DESIGN FINAL DESIGN 



NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

170.0072 

196.3611 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

30.0000 

0.0820 

3 

60.0000 

5.9920 

D 

90.0000 

89.9130 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

STRESS 
FAILURE 
CRITERIA 
F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F ■ 1 
FALURE : F > 1 

MATERIAL SYSTEM § 

1 

0.0172 

0.0183 

2 

0.0188 

0.0183 

3 

0.0221 

0.0184 

B 

0.0238 

0.0257 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

MAX. STATIC DEFLEC. 
IN 'X' DIREC. (IN.) 

3.50E-4 <A < 5.50E-4 

5.1486E-4 

3.5000E-4 
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CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.2 : TEST CASE #2 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

|— X 

y 1 

800 LBS. ' 

1 

T300 / IMHS 
FVR: 0.62 

TEMP. EACH NODE £ \ 

T U =250.0°F * 

OPTIMIZATION VARIABLES 

INITIAL DESIGN FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

170.0072 

194.4557 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

30.0000 

0.0820 

3 

60.0000 

16.7150 

H 

90.0000 

89.8920 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

STRESS 

2 

1 

0.0172 

0.0146 

2 

0.0188 

0.0146 


FAILURE 

LU 

*— 

GO 

3 

0.0221 

0.0199 

CRITERIA 

&> 
\ 

B 

0.0238 

0.0307 

| r * i.u- 

MUUL 

< 




NO FALURE : F < 1 

c r 

UJ 

5 

NOT USED 

NOT USED 

INUlrlLm 

FALURE 

: r « r 

:F > 1 

< 

2 

6 

NOT USED 

NOT USED 




B 

NOT USED 

NOT USED 


MAX. STATIC DEFLEC. 
IN y DIREC. (IN.) 

0.90E-2 <A < 1.15E-2 


0.9747E — 2 


1. 1500E-2 





























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.3 : TEST CASE #3 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 


T300 / IMHS 
FVR: 0.62 


OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

170.0072 

193.2950 

UtblUN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

30.0000 

0.0820 

3 

60.0000 

21.0780 

D 

90.0000 

89.9390 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

CUNbl KAINI 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDEIE 

NO FALURE : F < 1 
INCIPIENT : F * 1 
FALURE : F > 1 

MATERIAL SYSTEM jjf 

B 

0.0172 

0.0566 

2 

0.0188 

0.0357 

3 

0.022 1 

0.0264 

B 

0.0238 

0.0371 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

MAX. STATIC DEFLEC. 
IN DIREC. (IN.) 

2.50E-1 <A £ 3.00E-1 

2.8657E— 1 

2.4932E- 1 

\ 
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TABLE 5.4 : TEST CASE #4 


FIBER /MATRIX SYSTEM 

LOAD CONDITION ,m°nr L ?f) 1 


T300 / IMHS 
FVR: 0.62 

— - q D 1) 1 

TEMP. EACH NODE £ g 

T = 250.0° F ^ Zi 100 LBS. 

1 x (NODE 55) 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 1 FINAL DESIGN 


OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

170.0072 

180.4136 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM ft 

1 

0.0000 

0.0000 

2 

30.0000 

15.7390 

3 

60.0000 

51.4950 

D 

90.0000 

90.0000 

5 

NOT USED 

NOT USED . 

6 

L 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDEIE 

NO FALURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM ft 

i 

0.0172 

0.0170 

2 

0.0188 

0.0193 

3 

0.0221 

0.0252 

B 

0.0238 

0.0259 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

MAX. STATIC DEFLEC. 
IN 'T DIREC. (IN.) 

2.00E-2 < A <2.70E-2 

2. 1826E-2 

2.7027E-2 





























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


t 


TABLE 5.5 : TEST CASE #5 


FIBER/MATRIX SYSTEM 

LOAD CONDITION 

NO MECHANICAL 
LOADS APPLIED 

T300 / IMHS 
FVR: 0.62 

TEMP. EACH NODE 
T U =250.0°F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

ObJLCI IVL 
FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

170.0072 

177.0385 

ULOIUIN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

30.0000 

20.7020 

3 

60.0000 

54.4550 

D 

90.0000 

90.0000 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

UUINDIRMIIN 1 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDEIE 

NO FALURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0172 

0.1641 

2 

0.0188 

0.0175 

3 

0.0221 

0.0221 

B 

0.0238 

0.0250 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

MAX. STATIC DEFLEC. 
IN 'X' DIREC. (IN.) 

0.30E-4 <A < 1.50E-4 

0.3979E— 4 

0.3000E— 4 
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A detailed inspection of the variable values in Table 5.1 (test 
case No. 1) will show the anticipated trend associated with the particular 
load condition shown. The fiber direction should align itself in the same 
direction as the bending modulus for a said bending modulus to be optimal. 
Therefore, if the modulus in question is "611", the ply angle of each 
material system should tend toward zero degrees off the 1-1 axis. Closer 
examination of the ply angles shows that some stopped short of that zero 
degree target. The reason for this anomaly comes clear by scrutinizing 
the constraint variables. The stress failure criteria of each material 
system were well within the given limits, yet the limiting value of the 
maximum static deflection was reached and the design obtained was 
considered to be optimal. 

Tables 5.2 and Table 5.3 (test cases No. 2 and No. 3 
respectively), are very similar to the first validation case discussed 
above. The obvious differences are the load conditions and the direction 
of the maximum allowable deflections. Test case No. 2 has a uniform 
mechanical load applied in the y-direction, while test case No. 3 has a 
uniform load applied in the z-di recti on. In each case, the final design 
of the composite structure (as in test case No. 1) was determined by the 
maximum allowable deflection in the direction of the loads. 

The final two validation cases. No. 4 and No. 5, shown in Tables 
5.4 and 5.5, reassured this researcher that given a twisting load as well 
as a pure thermal load applied to the structure, the program will perform 
as expected. These two test cases revealed little when examined alone, 
but in concert with the other validation cases, they proved to be 
informative. 
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The initial design of all five validation cases was the same, but 
the final designs were determined by their individual load conditions and 
applied constraints. Since the same objective function is being optimized 
for each case, a comparison of the final design strength to ply angle 
directions can be made. It can be observed that the increase in the "Gil" 
NASTRAN bending modulus of the composite structure is in direct proportion 
to how close the ply angles came to the 1-1 axis. Though this can be 
regarded as obvious to the trained eye, it is reassuring to see this basic 
trend come to light. 

From the results of the test cases shown on Table 5.1 through 5.5, 
the modules contained in this program can be considered to have been 
successfully coupled. The depth of the validation outlined above was 
determined by the needs of the demonstration cases. The next chapter will 
expound on further recommendations for program expansion and increased 
validation testing. 

5.2 Demonstration Cases 

In Chapter 2, the fiber/matrix systems used for the demonstration 
cases were described in detail. This section will elaborate on the load 
conditions applied to each composite system chosen. The effects these 
different mechanical and thermal loads will have on the composite 
structure will determine the next step in establishing an optimal design 
for the aforementioned loading conditions. 

Tables 5.6 through 5.35, found in the Appendix (pages 72 to 106), 
contain the details of each demonstration case and will be referred to 
periodically throughout the rest of the chapter. In the title of these 
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tables, the demonstration case number also refers to the load condition 
number shown in figure 5.2; the letter adjacent to the load condition 
number refers to the following fiber/matrix systems: 

A: AS/IMHS FVR -0.60 

B: HMSF/IMHS FVR « 0.60 

C: HMSF/IMHS//SGLA/IMHS FVR -0.60 (each) 

For example, demonstration case "4C" will use load case "4" applied to 
fiber/matrix system "C". For quick reference, each table will contain the 
fiber/matrix system used as well as the load condition applied. 

As mentioned above, figure 5.2 shows the load conditions applied 
to each composite structure. The details of the mechanical and thermal 
loads are conveyed in a concise manner in this figure. A steady-state 
temperature of 250*F will be used with all six mechanical loading 
conditions. This steady-state condition is needed to isolate the effects 
of the variety of applied mechanical loads on the structure. The balance 
of the loads shown in figure 5.2 will be the four thermal gradients 
applied to the structure. These thermal conditions were designed to 
segregate thermal gradients in the three directions of the main axes as 
well as a gradient applied simultaneously in the same three directions. 

As stated in the previous section, the only constraints of concern 
for the demonstration cases are the stress failure criteria of each ply. 
This has been determined to be an adequate test as to the ability of the 
structure to withstand the load imposed upon it. If the ply stress 
failure criteria is less than 1.0, there is no ply failure, and therefore 
no failure of the structure. If the failure criteria is exactly 1.0, the 
ply is on the verge of failing, and when the stress failure criteria goes 
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OBJECTIVE 

FUNCTION 


x HZ 800 LBS. 

G1 1 (ELEMENT 12/13) 


FIXED 

TEMP. AT 
THIS FACE 
T U =300°F 



_ 7 _ 


TEMP. AT 
THIS FACE 
T U =70°F 


funct»on : G1 1 (ELEMENT 12/13) 


TEMP. AT THIS FACE Q 
T = 1 00 0 F ® 

^iv/r*r\ 1 (j — 


FIXED 

Zi 


OBJECTIVE 

FUNCTION 


f 800 LBS. 
G1 1 (ELEMENT 9/16) 


y 



X 


La 


L x T u =200°F 

TEMP. AT THIS FACE 

SEE : G11 (ELEMENT 12/13) 


TEMP. AT THIS FACE T u = 300 0 F 

\ yj 

FIXED I 


OBJECTIVE 

FUNCTION 


\ 800 LBS. 
G12 (ELEMENT 12/13) 


Z | x i 800 LBS. 


_ 4 _ 


OBJECTIVE 

FUNCTION 


G1 1 (ELEMENT 12/13) 


Z, 100 LBS. 

1 y i (NODE 51) 


100 LBS. 
(NODE 55) 


OBJECTIVE 

FUNCTION 


G1 1 (ELEMENT 9/16) 

Zi .100 LBS. A 

I x i (NODE 51) Lr_ 


FIXED I 

\ i I 

TEMP. AT THIS FACE T u = 100° F 
SKSEf: G11 (ELEMENT 12/13) 


T u = 70 ° F (MODE 11°) 10 

r~y i 

FIXED I 

A hn 

T u = 300° F (NODE 1) 

?Sic C T SJ : G11 (ELEMENT 12/13) 


NOTES: 

• TEMP. OF EACH NODE FOR 
LOAD CASES #1 THRU #6 : 

T u = 250.0° F 

• TEMP. OF EACH NODE FOR 
LOAD CASES #7 THRU #10 
ARE AS SHOWN 


100 LBS. 
(NODE 55) 


OBJECTIVE 

FUNCTION 


G12 (ELEMENT 12/13) 


Figure 5.2: Load Conditions 
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above 1.0, the ply has failed. At the point where an individual ply has 
failed, the entire structure is considered to have failed. 

The balance of this section will discuss some of the more 
interesting demonstration cases found on Tables 5.6 through 5.35. The 
ones chosen are considered to have problems that warrant discussion and 
the unique results are due to their particular combination of load 
condition, fiber/matrix system and optimization variables. 

A detailed look at Table 5.10 reveals the failure of material 
system No. 4 within the HMSF/IHMS composite structure that has load 
condition No. 2 imposed on it. Examining the ply angle of said material 
system shows that it is well off an expected play angle to given an 
optimal H G11" NASTRAN bending modulus. A second case. Table 5.10.1, was 
run using the same composite system and load condition, but with different 
design variables. In this case, the thickness of each material system, 
rather than the ply angle, was allowed to vary. 

Examining the results of this case shows that the optimal bending 
modulus was obtained by the proper distribution of the material system 
thicknesses to the system with the critical ply angle. Comparing the 
results of the two cases shown in Tables 5.10 and 5.10.1 confirms that the 
original bending modulus had reached an optimum value, within a very small 
percentage, the failure criteria however, had improved to a value well 
within the limits of failure. This demonstration case is proof of the 
sensitivity of the failure criteria to ply angle direction. 

The demonstration case that will be addressed next is shown on 
Table 5.16. The HMSF/IMHS composite structure failed with the initial 
design having load condition No. 4 bearing on it. It was observed, using 


61 



engineering judgement, that the structure failed due to inadequate overall 
composite thickness. This posed an interesting problem of finding the 
optimal thickness of the composite and still obtaining a maximum value for 
the M G11 M bending modulus. Table 5.16.1 shows the results of how this 
particular problem was attempted to be solved. The objective function 
used is the summation of the bending modulus and the Inverse of the 
composite thickness. Analysis of the results discloses an attempt by the 
program to vary the ply angles to increase the bending modulus, and at the 
same time, decrease the thickness to an optimal value. Before any 
effective change of the bending modulus could be made by changing the ply 
angle direction, the composite thickness was reduced to where the 
outermost plies were failing under the load. This is evident by the 
stress failure criteria of material system No. 1 increasing to a value 
above incipient failure and thus stopping the design. 

Using the final design thickness from Table 5.16.1, a third run 
was made with the same load conditions. This case used the objective 
function and design variables of the original case shown on Table 5.16. 
The outcome of this endeavor is displayed in Table 5.16.2. It can be 
seen that using the new composite thickness, a final optimal design was 
achieved which kept the stress failure criteria well below the value 
needed to fail the structure. 

The problem encountered using the HMSF/IMHS composite structure 
was repeated when the same load condition No. 4 was applied to the 
HMSF/ IMHS/SGLA/ IMHS intraply composite system. The same methodology as 
above was applied to this case top obtain a desired design that would not 
fail with the load bearing on it. The effect of using the problem solving 
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techniques formulated with the previous case is shown in greater detail 
on Tables 5.17, 5.17.1 and 5.17.2. The aspects of the final design are 
much the same as those shown in the preceding problem, yet with one big 
difference; the optimal thickness for this fiber/matrix system is smaller. 

5.3 S mnwiry 

Each load case was run with identical initial design variables and 
allowed to reach an optimal design dependent on the respective objective 
functions and constraints. The ply stress failure values were output for 
each case, and the maximum value is displayed in figure 5.3. The 
importance of this constraint variable follows from the fact that the 
failure of the entire structure hinges on its value. When addressing 
these results, the structural capacities inherent to each fiber/matrix 
system are revealed. This is evident from the effects of different load 
conditions on each composite system. The intrinsic capabilities of these 
composites can be shown to be directly attributed to the properties of the 
fibers used in the structure. This observation is based on the knowledge 
that the same matrix material and fiber volume ratio were used in each 
composite. Therefore, the differences encountered between the composite 
system are based strictly on the differences in fiber properties. 

Further reinforcement of this statement can be found by a closer 
examination of figure 5.3. The HMSF/IMHS composite structure appears to 
perform poorly compared to the other fiber/matrix systems. The maximum 
stress failure value of the HMSF/IMHS system is greater than the other two 
systems, yet the NASTRAN bending modulus is the largest. This apparent 
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lack of strength is traced to the fiber properties shown in figure 5.4. 
The strength of the "AS" and "SGLA" fibers is larger than the "HMSF", but 
the normal moduli is lower. Taking this into consideration, it is now 
understood why the overall composite strength was improved when the "SGLA" 
fiber was intermixed with the "HMSF" fiber. 


65 



Name : Symbol 

Units 

AS 

SGLA 

HMSF 

Density 

f 

lb /in 3 

0.063 

0.090 

0.070 

Longit. modulus 

E 

f 11 

10 6 psi 

31.00 

12.40 

55.00 

Transv. modulus 

E 

f22 

10 6 psi 

2.000 

12.40 

0.900 

Long, shear mod. 

G 

f 12 

10 6 psi 

2.000 

5.170 

1.100 

Trans, shear mod. 

: G 

f23 

10 6 psi 

1.000 

5.170 

0.700 

Longitudinal 
Poisson's ratio 

f 12 

— 

0.200 

0.200 

0.200 

Transverse 
Poisson's ratio 

f23 

— 

0.250 

0.200 

0.250 

Heat capacity 

C 

f 

Btu 

0.200 

0. 170 

0.200 

■na 

Longitudinal 
heat conductivity 

: K 

f 1 1 

IHI 

580.0 

21.00 

i 

580.0 

hr-ft 2 -°F-in 

Transverse 

:K 

f22 

Btu 

58.00 

i 

21.00 

58.00 

heat conductivity 

hr- ft 2 -°F-in 

Longit. thermal 
expansion coeff. 


in" 6 in 

-55.0 

2.800 

-55.0 

f 1 1 

1 U . o 

in - F 

Transv. thermal 
expansion coeff. 


in -6 in 

5.600 

2.800 

5.600 

f22 

10 . Or- 

in - F 

Longitudinal 
tension strength 

:S ft 

ksi 

350.0 

600.0 

250.0 

Longitudinal 
compressive str. 

:S fc 

ksi 

260.0 

540.0 

200.0 


Figure 5. 4 : Fiber Properties 
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CHAPTER VI 


OBSERVATIONS, RECOMMENDATIONS AND SUMMARY 

As researchers strive to improve their understanding of new and 
existing materials, the marriage of laboratory testing and computer 
software development was a necessity. From this union, many different 
disciplines of study were able to develop accurate and capable programs. 
As computer software development evolved, it was evident that a more "real 
life" loads environment could be duplicated without the major expense of 
repeated labor and hardware intensive laboratory test. I am not 
advocating abandoning the use of laboratory testing, but with the 
application of multi discipline software, expensive testing processes can 
be reduced. This can be accomplished by generating a design from a 
verified multidiscipline program. This design will be tailored to a 
combined load condition derived from surroundings the structure is exposed 
to. The final step is to confirm the design by duplicating the 
environment in a test cell and testing a prototype of the structure. 

This project could be considered to be a step toward the evolution 
of multi discipline software in the field of composite mechanics. With the 
successful integration of programs involving composite mechanics, heat 
transfer, structural analysis and structural optimization, it becomes 
evident that a composite structure can be designed to a combined loads 
environment. This first step individually applied various thermal and 
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mechanical loads to a composite structure. The computer code notably used 
ply angles, ply thickness and combined stress criteria to affect the 
design of the composite. The integrated program performed well in 
designing a composite upon which different isolated mechanical loads were 
imposed. The effects of combined mechanical load conditions were not 
detailed in this research. Combinations of the load cases used in this 
research are left to be run as a follow-on to this research. The 
evaluation of results from the combined mechanical load states will lead 
to a more complete study of the fiber/matrix systems used. 

Elaborating on this statement, it is evident that the expansion 
of this research is necessary. As mentioned above, increasing the variety 
of demonstration cases is one suggestion. Incorporating different 
composite properties as design variables could greatly increase the number 
of design possibilities. The use of cost or weight constraints could be 
helpful since these are two important design drivers in the space program. 
Examining the other fields of study, such as the structural or thermal 
disciplines, would also expand the possibilities for more research. The 
structures area could include dynamic loading, not just the static loading 
investigated here. In the region of heat transfer, the thermal distortion 
of a structure could be taken into consideration. Fatigue life cycles of 
a composite structure can also be studied. 

In all examples cited above, a full spectrum of validation testing 
will be needed. As each milestone is reached, verification testing should 
follow as a means of independent demonstration that true results were 
obtained. 
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To reinforce the conclusions, and since the main thrust of this 
research was the effects of composite materials in a multidiscipline 
environment, any further research should start here. Expansion of the 
optimization variables to include the fiber and matrix properties used in 
ICAN, and including the fiber volume ratio, could help if weight and/or 
cost versus strength became an issue. 
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CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.6 : DEMO. CASE # 1A 


fTBER/MATRIX SYSTEM 
AS 7 IMHS 
FVR: 0.60 


TEMP. EACH NODE 


T u = 250.0 F 


OPTIMIZATION VARIABLES INITIAL DESIGN 



FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

175.0035 

188.0425 



i 

0.0000 

0.0000 



2 

18.0000 

- 0.1300 

PLY 

LJ 

h- 

00 

3 

36.0000 

0.4320 

ANGLES 

>- 

00 

-J 

H 

54.0000 

‘ 46.2840 

(DEGREES) 

< 

cr 

LJ 

5 

NOT USED 

NOT USED • 


5 

6 

NOT USED 

NOT USED 



7 

NOT USED 

NOT USED 



1 

0.0043 

0.0023 

STRESS 

FAILURE 


2 

0.01 13 

0.0024 

LJ 

1— 

CO 

3 

0.0253 

0,0026 

CRITERIA 

>- 

CO 

B 

0 0375 

0 0453 

BBQ3BS33 

1 


V.VT JJ 

< 




NO FAILURE : F < 1 
INCIPIENT : F « 1 
FAILURE : F > 1 

cr 

LJ 

5 

NOT USED 

NOT USED 

!<., 

2 

6 

NOT USED 

NOT USED 



B 

NOT USED 

NOT USED 
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CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.7 : DEMO. CASE # IB 


FIBER/ ^.MATRIX SYSTEM .LqAD.CONpmON Q J 
HMSF/IMHS TEMP. EACH NODE xj 

FVR: 0.60 T =250.0°F u '' 


T u = 250.0 F 


OPTIMIZATION VARIABLES INITIAL DESIGN 


FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 



=«= 2 

2 

uj , 
ti 3 


cr 5 
uj ° 


=*| 2 
UJ 


STRESS 

FAILURE 

CRITERIA 


NO FAILURE : F < 1 5 

INCIPIENT : F - 1 H- 
FAIURE : F > 1 S' 


308.0044 

331.8291 

0.0000 

0.0000 

18.0000 

-0.2684 

36.0000 

0.3800 

54.0000 

38.4180 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

0.0035 

0.0996 

0.0586 

0.1055 

0.1371 

0.1080 

0.1553 

0.8154 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 
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TABLE 5.8 : DEMO. CASE # 1C 


FIBER /MATRIX SYSTEM 

LOAD CONDITION p 

i 

X 

J 

>* 


HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

TEMP. EACH NODE i 
T u = 250.0° F u 

' ^ QJ 

-♦o 

L^o 

00 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

260.7607 

280.7705 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

2 

18.0000 

0.0610 

3 

36.0000 

0.5410 

H 

54.0000 

19.0430 

5 

NOT USED 

NOT USED ' 

6 

NOT USED 

NOT USED 

7 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FAILURE : F < 1 
INCIPIENT :F ■ 1 
FALURE : F > 1 

MATERIAL SYSTEM § 

n 

0.0049 

0.0007 

2 

0.0199 

0.0007 

3 

0.0448 

0.0010 

D 

0.0662 

0.0395 

5 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

























































TABLE 5.9 : DEMO. CASE #2A 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

o i 

Ld 1 
* 

y | 

X 

800 LBS. i 


AS / IMHS 
FVR: 0.60 

TEMP. EACH NODE 
T u = 250.0° F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 9/16) 

175.0035 

188.0425 

i 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM § 

1 

0.0000 

0.0000 

2 

18.0000 

- 0.1300 

ZLd 

3 

36.0000 

0.4320 

— CD 

</>< 

D 

54.0000 

46.2840 

5 

NOT USED 

NOT USED ' 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F * 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0043 

0.0103 

2 

0.01 13 

0.0104 

3 

0.0253 

0.0107 

B 

0.0375 

0.0519 

5 

NOT USED 

NOT USED 

6 

i 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 























































TABLE 5.10 : DEMO. CASE#2B 


FIBER/MATRIX SYSTEM 

LOAD CONDITION 

o 

LJ < 
* 

y l 

X 

800 LBS. i 


HMSF / IMHS 
FVR: 0.60 

TEMP. EACH NODE 
T u = 250.0° F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 9/16) 

308.0044 

331.6025 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM ff 

n 

0.0000 

0.0000 

2 

18.0000 

2.0560 

3 

36.0000 

0.9800 

D 

54.0000 

54.0000 

5 

NOT USED 

NOT USED • 

6 

NOT USED 

NOT USED 

H 

NOT USED 

NOT USED 

CONSTRAINT 
| VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 

MATERIAL SYSTEM ff 

n 

0.0035 

0.5153 

2 

0.0586 

0.5931 

3 

0.1371 

0.5510 

II 

0.1553 

1.3752 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 
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TABLE 5.10.1 : DEMO. CASE #2B2 


.QBER/ MAJRJ.X SYSTEM I LOAD CONDITION 
HMSF / IMHS 
FVR: 0.60 


OPTIMIZATION VARIABLES INITIAL DESIGN 


FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G11 
g| (ELEMENT 9/16) 


308.0044 


0.0500 © 00° 


=*= 2 0.0500 @ 18° 


MATERIAL 5 

CO ►— 

gy THICKNESS > 
^cq < J ) 

(INCHES ® ^ 


»< — 
> PLY ANGLE) £ 


K 5 


< 

2 6 


0.0500 @ 36° 


0.0050 @54° 


NOT USED 


NOT USED 


331.841 1 


0.1550 © 00° 


0.0000 © 18° 


0.0000 @36° 


0.0000 @54° 


NOT USED 


NOT USED 


NOT USED 


NOT USED 


STRESS 

FAILURE 

CRITERIA 


=*fc 2 


NO FAILURE 
INCIPIENT 
FALURE 


f <' 1 5 

f»’ IT 


0.0035 


0.0586 


0.1371 


0.1553 


NOT USED 


NOT USED 


0.2000 



NOT USED 


NOT USED 


NOT USED 




















































TABLE 5. 11 : DEMO. CASE # 2C 


% PLY: 80 % PLY: 20 


OPTIMIZATION VARIABLES 


UJ 2 


LOAD CONDITION 

o \ 


TEMP. EACH NODE 

Id 

2 

y |_ 

T U =250.0°F 

li. / 

A 

X 

800 LBS. < 


INITIAL DESIGN 

FINAL DESIGN 

260.7607 

280.7705 

0.0000 

0.0000 

18.0000 

0.0610 

36.0000 

0.5410 

54.0000 

19.0430 

NOT USED 

NOT USED ’ 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

0.0049 

0.0332 

0.0199 

0.0332 

0.0448 

0.0334 

0.0662 

0.0347 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 


£o NASTRAN BENDING 
£u MODULUS: G 1 1 
(ELEMENT 9/16) 


•M 2 

2 1 
LxJ 


CO 

ANGLES oo 


(DEGREES) ^ 


8 5 
< 

2 6 


STRESS 

FAILURE 

CRITERIA 


=**= 2 

2 

Ld - 

H- 3 


NO FALURE 
INCIPIENT 
FALURE 


F < 1 

LJ 

F-1 £ 






















































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.12 : DEMO. CASE #3A 


fJBER/MATRIX SYSTEM 
AS / IMHS 
FVR: 0.60 



OPTIMIZATION VARIABLES INITIAL DESIGN FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G12 
(ELEMENT 12/13) 

8.2800 

16.0392 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

2 

18.0000 

44.9737 

3 

36.0000 

46.1552 

D 

54.0000 

53.9851 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

7 

NOT USED 

NOT USED 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDEIE 

NO FAILURE : F < 1 
INCIPIENT : F ■ 1 
FAILURE : F > 1 

MATERIAL SYSTEM § 

n 

0.0043 

0.0133 

2 

0.01 13 

0.0249 

3 

0.0253 

0.0253 

B 

0.0375 

0.0209 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 


79 

























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


DEMO. CASE # 3B 


TABLE 5.13 : 


.EIBER/MAJRIX SYSTEM 
HMSF / IMHS 
FVR: 0.60 


OPTIMIZATION VARIABLES INITIAL DESIGN 



NASTRAN BENDING 
MODULUS: G12 
(ELEMENT 12/13) 


2 


12.2744 


0.0000 


18.0000 


FINAL DESIGN 


26.4908 


0.0000 


44.8047 


to 

ANGLES 


(DEGREES) | 

Ld ^ 
<: 

2 6 


=*=1 2 


STRESS 

FAILURE 

CRITERIA 


NO FALURE : F < 1 Ct 5 

INCIPIENT : F * 1 

FALURE : F > 1 < - g 


Dl 


36.0000 


54.0000 


NOT USED 


NOT USED 


NOT USED 


0.0035 


0.0586 


0.1371 


0. 1553 


NOT USED 


NOT USED 


NOT USED 


45.9202 


53.9840 


NOT USED 


NOT USED 


NOT USED 


0.0072 


0.0288 


0.0287 


0.0256 


NOT USED 


NOT USED 


NOT USED 
























































TABLE 5.14 : DEMO. CASE#3C 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

o \ 

ui 

>< 

y l 

X 

800 LBS. . 


HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

TEMP. EACH NODE 
T u = 250. 0°F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G12 
(ELEMENT 12/13) 

10.7768 

22.6956 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

45.0151 

3 

36.0000 

46.1603 

D 

54.0000 

53.9770 

5 

NOT USED 

NOT USED • 


6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0049 

0.0165 

2 

0.0199 

0.0356 

3 

0.0448 

0.0360 

D 

0.0662 

0.0283 

5 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 






















































TABLE 5.15 : DEMO. CASE #4A 


FIBER/MATRIX SYSTEM 

LOAD CONDITION 

Q 800 LBS. k 

X 3 

Z |— X 

AS / IMHS 
FVR: 0.60 

TEMP. EACH NODE 
T U =250.0°F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 11 
(ELEMENT 12/13) 

175.0035 

188.0326 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

- 0.3961 

3 

36.0000 

1.1737 

B 

54.0000 

50.1540 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F- 1.0- MDEIE 

NO FALURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0043 

0.4846 

2 

0.01 13 

0.2169 

3 

0.0253 

0.0568 

B 

0.0375 

0.0516 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 


82 



























































TABLE 5.16 : DEMO. CASE #4B 


fiBER/MATRIX system load condition q 

HMSF / IMHS TEMP. EACH NODE H 

FVR: 0.60 


800 LBS. 


T u = 250.0 F 


OPTIMIZATION VARIABLES INITIAL DESIGN 



FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

308.0044 



B 

0.0000 



2 

18.0000 

PLY 

UJ 

CO 

3 

36.0000 

ANGLES 

CO 

_J 

B 

54.0000 

(DEGREES) 

< 

cc 

UJ 

5 

NOT USED 



6 

NOT USED 



B 

NOT USED 



1 

8.0508 

QTDpCC 


2 

3.1547 

! O l r\ O O 

FAILURE 

UJ 

IT) 

3 

0.6581 

CRITERIA 

fc 

B 

0 1687 


1 


£ 



NO FAILURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

cr 

UJ 

5 

NOT USED 

<. 

2 

6 

NOT USED 



B 

NOT USED 



NOT USED 


NOT USED 


NOT USED 



NOT USED 


NOT USED 


NOT USED 


83 













































TABLE 5. 16.1 : . DEMO. CASE #4B2 


FIBER/ MATRIX SYSTEM 

LOAD CONDITION 

Q 800 LBS. i 

Vi— - 1 

Ll. 

Z L X 

HMSF / IMHS 
FVR: 0.60 

TEMP. EACH NODE 
T u = 250.0° F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

BENDING INVERSE OF 

MODULUS 4 " THICKNESS 
( G 1 1 + 1.0 / t c ) 

(ELEMENT 12/13) 

244.2977 

261.7686 


PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

DESIGN 

VARIABLES 

2 

18.0000 

17.6872 

3 

36.0000 

35.0955 

fl 

54.0000 

53.4269 

5 

72.0000 

71.9281 • 

6 

90.0000 

90.0000 

E 

0.0000 

l 

l 

1 

l 

COMPOSITE / IM \ 
THICKNESS 

0.6200 

0.5484 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT :F-1 
FAILURE : F > 1 

MATERIAL SYSTEM # 

i 

0.01 16 

1.0015 

2 

0.0117 

0.5486 


0.0118 

0.1978 

D 

0.0120 

0.0644 

5 

0.0122 

0.0254 

6 

0.0122 

0.0148 

E 

0.0116 

1 

1 

1 

1 


84 























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.16.2 : DEMO. CASE #4B3 


FIBER /MATRIX SYSTEM 

LOAD CONDITION Q 800 LBS. | 

HMSF / IMHS 
FVR: 0.60 

TEMP. EACH NODE x § 

T =250.0° F ^ z. 

1 X 

OPTIMIZATION VARIABLES 

INITIAL DESIGN FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

257.1384 

331.5830 



n 

0.0000 

0.0000 



2 

18.0000 

0.1660 

PLY 

LJ 

H- 

(/) 

3 

36.0000 

0.1628 

ANGLES 

>“ 

00 

_l 

D 

54.0000 

0.4220 

(DEGREES) 

< 

oc 

LJ 

5 

72.0000 

1.3700 


< 

6 

90.0000 

95.0850 



7 

NOT USED 

NOT USED 



n 

0.0110 

0.5163 

STRESS 


2 

0.01 13 

0.3287 

FAILURE 

LJ 

♦— 

00 

3 

0.01 19 

0.1836 

CRITERIA 

> 

00 

B 


0.0808 

E3ESE233 

1 

0.0 125 

< 




NO FALURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 

cr 

LJ 

5 

0.0128 

0.0208 

i 

6 

0.0130 

0.0476 



B 

NOT USED 

NOT USED 


85 

























































TABLE 5.17 : DEMO. CASE #4C 


fJBER/MATRIX SYSTEM L OAD CONDiHON Q 

hmsf/imhs^sgla/imhs temp. EACH NODE xC 

FVR: 0.60 FVR: 0.60 j = 250 0°F u " 

% PLY: 80 % PLY: 20 u 


OPTIMIZATION VARIABLES INITIAL DESIGN 


800 LBS. 


NASTRAN BENDING 
MODULUS: G11 
§| (ELEMENT 12/13) 


2 


260.7607 


0.0000 


18.0000 


in 

ANGLES in 


(DEGREES) | 

i . i O 


36.0000 


54.0000 


NOT USED 


NOT USED 


NOT USED 


3.4189 


STRESS 

FAILURE 

CRITERIA 


NO FALURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 


=*| 2 


1.1535 


0.1059 


0.0678 


NOT USED 


NOT USED 


NOT USED 


NOT USED 


NOT USED 



86 















































TABLE 5.17.1 : DEMO. CASE #4C2 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

Q 800 LBS. i 

>< 3 

Z L* 


HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

TEMP. EACH NODE 
T U =250.0°F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

BENDING INVERSE OF 

MODULUS + THICKNESS 
( G 1 1 + 1.0 / t c ) 

(ELEMENT 12/13) 

207.5044 

247.7756 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

— - 

1 

0.0000 

0.0000 

2 

18.0000 

16.7065 

3 

36.0000 

32.6863 

D 

54.0000 

52.2459 

5 

72.0000 

71.8550 

6 

90.0000 

— 


D 

0.0000 

— 

COMPOSITE / 1 ki \ 
THICKNESS u ' 

0.6200 

0.4270 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F~ 1.0 — MDEIE 

NO FALURE : F < 1 
INCIPIENT : F ■ 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0262 

0.9977 

2 

0.0263 

0.4356 

3 

0.0266 

0.0849 

D 

0.0269 

0.0439 

5 

0.0271 

0.0423 

6 

0.0272 

— 

D 

0.0262 

— 






















































TABLE 5.17.2 : DEMO. CASE #4C3 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

Q 800 LBS. i 

x t 

Z L- x 

HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

TEMP. EACH NODE 
T U =250.0°F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

240.9406 

280.7046 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

0.2610 

. .. . j 

3 

36.0000 

- 0.4880 

H 

54.0000 

1.0140 

5 

72.0000 

71.6960 • 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

I CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FAILURE : F < 1 
INCIPIENT : F * 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

n 

0.0180 

0.6513 

2 

0.0216 

0.3578 

3 

0.0267 

0.1514 

B 

0.0328 

0.0350 

5 

0.0376 

0.0980 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 


88 
























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.18 : DEMO. CASE # 5A 


.E'BER/MATRIX SYSTEM | LOAD .CONDITION 
AS/IMHS 
FVR: 0.60 


100 LBS. 
(NODE 51) 


OPTIMIZATION VARIABLES INITIAL DESIGN 


FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 9/16) 


=*= 2 
Ld 


CO 

ANGLES in 


(DEGREES) g 

uj 

*< 


cn 5 
uj ° 


STRESS 

FAILURE 

CRITERIA 


=*| 2 

2 
UJ 


NO FAILURE 
INCIPIENT 
FALURE 


£ 

F>1 i 6 


175.0035 

188.0326 

0.0000 

0.0000 

18.0000 

- 0.1300 

36.0000 

0.4320 

54.0000 

46.2840 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

0.0043 

0.0314 

0.01 13 

0.0182 

0.0253 

0.0081 

0.0375 

0.0466 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 


89 






















































TABLE 5.19 : DEMO. CASE#5B 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

100 LBS., 
O (NODE 5 1) 1 

Ld a 

*3 

^ Z| 100 LBS. 

1 x (NODE 55) 


HMSF / IMHS 
FVR: 0.60 

TEMP. EACH NODE 
T U =250.0°F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 9/16) 

308.0044 

331.4558 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

2 

18.0000 

- 1.8559 

3 

36.0000 

4.74 15 

H 

54.0000 

53.1936 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

H 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDElE 

NO FALURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0035 

1.0014 

2 

0.0586 

0.7041 

3 

0.1371 

0.4648 

B 

0.1553 

0.8008 

5 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 


90 






















































TABLE 5.20 : DEMO. CASE#5C 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

100 LBS. i 
q (NODE 51)' 

UJ 

X 3 

^ Z| 100 LBS. < 

1 x (NODE 55) 


HMSF/IMHS/^SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

TEMP. EACH NODE 
T u = 250.0° F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 9/16) 

260.7607 

280.7705 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

; 

0.0000 

0.0000 

2 

18.0000 

0.0610 

i 

3 

L 

36.0000 

0.5410 

D 

54.0000 

19.0430 

5 

NOT USED 

NOT USED • 

6 

i 

NOT USED 

NOT USED 

H 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDEIE 

NO FALURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

n 

0.0049 

0.0588 

2 

0.0199 

0.0289 

3 

0.0448 

0.0097 

n 

0.0662 

0.0416 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 


91 























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.21 : DEMO. CASE # 6A 


f.'BlR/MATRIX .SYSTEM 
AS / IMHS 
FVR: 0.60 


LOAD CONDI TION 
TEMP. EACH NODE 
T u = 250.0° F 


100 LBS. 

O (NODE 5 1) 1 

LlJ * 

* 3 

^ 2| 100 LBS. i 

I X (NODE 55) 


OPTIMIZATION VARIABLES INITIAL DESIGN 


FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G12 
(ELEMENT 12/13) 

8.2800 

16.0392 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

44.9737 

3 

36.0000 

46.1552 

D 

54.0000 

53.9851 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

STRESS 

FAILURE 

CRITERIA 

F« 1.0-MDEIE 

NO FALURE : F < 1 
INCIPIENT : F - 1 
FAILURE : F > 1 

MATERIAL SYSTEM § 

n 

0.0043 

0.0276 

2 

0.0113 

0.0332 

3 

0.0253 

0.0270 

II 

0.0375 

0.0234 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 


























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.22 : DEMO. CASE # 6B 


FIBER /MATRIX SYSTEM 

LOAD CONDITION '“ LBS. < 


HMSF / IMHS 
FVR: 0.60 

jj ^INUUL j I j 

TEMP. EACH NODE " 

T U =250.0 U F ^ 2i 100 LBS.- 

1 x (NODE 55) 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 1 FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G12 
(ELEMENT 12/13) 

12.2744 

26.4905 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM jf 

1 

0.0000 

0.0000 

2 

18.0000 

44.8097 

3 

36.0000 

45.9582 

D 

54.0000 

53.9830 

5 

NOT USED 

NOT USED - 

6 

NOT USED 

NOT USED 

.7 

NOT USED 

NOT USED 

STRESS 

FAILURE 

CRITERIA 

F« 1.0-MDEIE 

NO FAILURE : F < 1 
INCIPIENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM § 

1 

0.0035 

0.0901 

2 

0.0586 

0.0937 

3 

0.1371 

0.0515 

D 

0.1553 

0.0253 

5 

NOT USED 

NOT USED 

m 

NOT USED 

NOT USED 





NOT USED 


NOT USED 





















































TABLE 5.23 : DEMO. CASE # 6C 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

100 LBS. J 

O (NODE 5 1) 

LJ s 

xf 

Z| 100 LBS. 

1 x (NODE 55) 


HMSF/ IMHS //SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

TEMP. EACH NODE 
T u « 250.0° F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G12 
(ELEMENT 12/13) 

10.7767 

22.6956 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM ff 

1 

I 

0.0000 

0.0000 

2 

18.0000 

45.0151 

3 

36.0000 

46.1603 

D 

54.0000 

53.9770 

5 

NOT USED 

NOT USED • 

6 

NOT USED 

NOT USED 

H 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 

MATERIAL SYSTEM ff 

1 

0.0049 

0.0471 

2 

0.0199 

0.0519 

3 

0.0448 

0.0400 

D 

0.0662 

0.0324 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 
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CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.24 : DEMO. CASE#7A 


TJB.ER/MAJRJ,X SXSTEM 
AS / IMHS 
FVR: 0.60 


LOAD CONDITION temp. 4 
THIS 4 



FIXED 


OPTIMIZATION VARIABLES 


NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 


PLY 

ANGLES 

(DEGREES) 


STRESS 

FAILURE 

CRITERIA 


F« 1.0- MDEIE 


NO FAILURE 
INCIPIENT 
FALURE 


F < 1 
F - 1 
F > 1 




LO 

>- 

CO 


< 

o: 




00 

> 

00 


< 

or 

LlJ 


1 


INITIAL DESIGN 

FINAL DESIGN 

174.9219 

187.9811 

0.0000 

0.0000 

18.0000 

- 0.1956 

36.0000 

- 0.4041 

54.0000 

30.9430 

NOT USED 

NOT USED • 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

0.0094 

0.0015 

0.0237 

0.0016 

0.0512 

0.0018 

0.0749 

0.0519 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 
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TABLE 5.25 : DEMO. CASE#7B 


HMSF / IMHS 
FVR: 0.60 


£o NASTRAN BENDING 
So MODULUS: G11 
S | (ELEMENT 12/13) 


=«= 2 

2 1 

LJ 


CO 

ANGLES £ 


(DEGREES) | 

UJ 5 
2 6 


STRESS 

FAILURE 

CRITERIA 


•fe 2 

2 
UJ 


NO FALURE 
INCIPIENT 
FALURE 


f < i ce 5 

F-1 
F > 1 < 


LOAD CONDITION TEMP. y temp. 

THIS ' y| THIS 

FACE / 1 x FACE 

3oo ° F kE 7o ° F 

INITIAL DESIGN 

FINAL DESIGN 

307.9368 

329.0864 

0.0000 

0.0000 

18.0000 

- 0.5280 

36.0000 

18.2170 

54.0000 

54.0000 

NOT USED 

NOT USED • 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

0.0073 

0.0399 

0.1149 

0.0464 

0.2665 

0.5062 

0.3015 

1.0000 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 
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TABLE 5.26 : DEMO. CASE #7C 


FIBER /MATRIX SYSTEM 
HMSF/IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

LOAD CONDITION temp. a 

.......................... TH1S j 

FACE > 
300 °F ' 

y L_ 

X 

FIXED 

TEMP. 

THIS 

FACE 

70 °F 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 


OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

260.6892 

280.7083 




1 

0.0000 

0.0000 



1 

2 I 
LJ 

00 

2 

18.0000 

- 0.1785 

00 

2LlJ 

PLY 

3 

36.0000 

0.0160 

o=! 

CD 

O0< 

ANGLES 

—J 

D 

54.0000 

54.0000 

(DEGREES) 

< 

cr 

LJ 

5 

NOT USED 

NOT USED . 



5 

6 

NOT USED 

NOT USED 




D 

NOT USED 

NOT USED 




n 

0.01 14 

0.0328 


STRESS 

i 

2 

0.0447 

0.0338 

zoo 

"Zf LJ 

FAILURE 

LJ 

CO 

3 

0.0953 

0.0329 


CRITERIA 

b\ 

D 

0.1387 

0.2614 






< 




o> 

NO FALURE : F < 1 
INCIPIENT : F ■ 1 
FALURE : F > 1 

Qd 

y ; 

5 

NOT USED 

NOT USED 


2. 

B 

NOT USED 

NOT USED 




B 

NOT USED 

NOT USED 
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TABLE 5.27 : DEMO. CASE#8A 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

TEMP. THIS FACE 
100° F 

1 J 

Z| 200° F 

1 — x temp. THIS FACE 

AS / IMHS 
FVR: 0.60 

FIXED 

Jill 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

175.5846 

188.4926 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

2 

18.0000 

0.0290 

3 

36.0000 

- 1.0090 

H 

54.0000 

54.0000 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F « 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F « 1 
FAILURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0079 

0.0163 

2 

0.0044 

0.0069 

3 

0.0037 

0.0015 

B 

0.0048 

0.0074 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 























































TABLE 5.28 : DEMO. CASE #8B 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

TEMP. THIS FACE 
100° F 

1 "ZD 

Zi 200° F 

1 x TEMP. THIS FACE 

HMSF / IMHS 
FVR: 0.60 

FIXEC 

ill! 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

308.4839 

332.2351 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

- 0.0470 

3 

36.0000 

0.3720 

B 

54.0000 

32.3690 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDElE 

NO FALURE : F < 1 
INCIPCNT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

B 

0.0537 

0.3944 

2 

0.0251 

0.1826 

3 

0.0361 

0.0482 

B 

0.0215 

0.0972 

5 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 





















































TABLE 5.29 : DEMO. CASE #8C 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

TEMP. THIS FACE 
100° F 

i 

Zi 200° F 

1 x TEMP. THIS FACE 

HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

FIXED 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

261.2646 

281.1792 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

. 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

2 

18.0000 

- 0.0194 

3 

36.0000 

0.5820 

D 

54.0000 

37.0820 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 
FAILURE 
CRITERIA 
F- 1.0- MDEIE 

MATERIAL SYSTEM § 

1 

0.0176 

0.0846 

2 

0.0043 

0.0388 

3 

0.0070 

0.0097 


0.008 1 

0.0141 

NO FALURE : F < 1 
INCIPIENT : F ■ 1 
FALURE : F > 1 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 


100 























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.30 : DEMO. CASE #9A 


FIBER/MATRIX SYSTEM 
AS / IMHS 
FVR: 0.60 

LOAD CONDITION 
o 

LJ 

X 

Li. 

300° F TEMP. THIS FACE 

y l_ 

i * ... ! 

100° F TEMP. THIS FACE 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

175.5099 

188.4361 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM jf 

n 

0.0000 

0.0000 

2 

18.0000 

- 0.0840 

3 

36.0000 

- 0.7670 

D 

54.0000 

51.0419 

5 

NOT USED 

NOT USED . 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 


STRESS 


n 

0.0076 

0.0040 

2 

0.0194 

0.0040 


FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 


UJ 

</) 

3 

0.0423 

0.0047 

>" 

00 

I 

B 

0.0620 

0.0850 

< 

cr 

u 

5 

NOT USED 

NOT USED 

!< 

2* 

6 

NOT USED 

NOT USED 


B 

NOT USED 

NOT USED 























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.31 : DEMO. CASE0 9B 


JJBER.Z MATRIX SYSTEM LOAD CO NDITION 300° f temp, t his face 

HMSF/IMHS y | 

FVR: 0.60 H *_ 

100° F TEMP. THIS FACE 


OPTIMIZATION VARIABLES INITIAL DESIGN FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 


=** 2 

2 ' 

Ld 


t n 

ANGLES 


(DEGREES) | 

Ld 

< 


on 5 

Ld 3 


STRESS 

FAILURE 

CRITERIA 


6 


7 


1 

2 


NO FALURE 
INCIPIENT 
FALURE 


3 

F < 1 OH 5 

F-1 £ — 
F> ’ 1 6 


308.422 1 

330.1665 

* 0.0000 

0.0000 

18.0000 

- 2.2610 

36.0000 

14.5310 

54.0000 

53.8850 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

0.0060 

0.0473 

0.0959 

0.0820 

0.2229 

0.4064 

0.2522 

0.9946 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 

NOT USED 


102 


























































TABLE 5.32 : DEMO. CASE#9C 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

300° F TEMP. THIS FACE 

T- 

x 

100° F TEMP. THIS FACE 

HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

Q 

UJ 

X 

u_ 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

261.1997 

281.1252 

</) 

2LlJ 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

n 

0.0000 

0.0000 

2 

18.0000 

0.2130 

3 

36.0000 

0.2600 

O— 1 

— QQ 
CO< 
LU~ 

D 

54.0000 

38.4200 

5 

NOT USED 

NOT USED • 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F« 1.0- MDElE 

NO FAILURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

1 

0.0091 

0.0101 

2 

0.0357 

0.0107 

3 

0.0775 

0.0109 

D 

0.1132 

0.1771 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

fl 

NOT USED 

NOT USED 


103 
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TABLE 5.33 : DEMO. CASE # 10A 


FIBER/MATRIX SYSTEM i 

iLOAD CONDITION 70° F TEMP.(NODE 110) 

AS / IMHS 
FVR: 0.60 

Q 

Ui 

X 

Ll. 

300° F TEMP, 


OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

NASTRAN BENDING 
MODULUS: G11 
(ELEMENT 12/13) 

175.4686 

188.3968 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

- 0.4310 

3 

36.0000 

- 0.5720 

D 

54.0000 

50.8310 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

E 

NOT USED 

NOT USED 


STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPENT : F - 1 
FALURE : F > 1 

MATERIAL SYSTEM # 

n 

0.0058 

0.0035 

2 

0.0143 

0.0036 

CONSTRAINT 

VARIABLES 

3 

0.0313 

0.0035 

B 

0.0460 

0.0553 


NOT USED 

NOT USED 

u 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 























































CONSTRAINT DESIGN OBJECTIVE 

VARIABLES VARIABLES FUNCTION 


TABLE 5.34 : DEMO. CASE # 10B 


FIBER/ MATRIX SYSTEM 

LOAD CONDITION 

IfSWSSmM 

MODE 110) 

HMSF / IMHS 
FVR: 0.60 

" o 

LjJ 

X 

u. 

300 °F TEMP.( 

OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 


NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

308.3884 

332.1479 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM # 

1 

0.0000 

0.0000 

2 

18.0000 

- 0.0780 

3 

36.0000 

0.5690 

D 

54.0000 

35.1000 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDEIE 

NO FALURE : F < 1 
INCIPENT : F - 1 
FALURE : F > 1 

1 

MATERIAL SYSTEM # 

B 

0.0045 

0.0837 

2 

0.0728 

0.0848 

3 

0.1686 

0.0951 

B 

0.1894 

0.8240 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 

























































TABLE 5.35 : DEMO. CASE # IOC 


FIBER /MATRIX SYSTEM 

LOAD CONDITION 

70 °F TEMP. (NODE 110) 

HMSF/ IMHS//SGLA/IMHS 
FVR: 0.60 FVR: 0.60 
% PLY: 80 % PLY: 20 

* — o 

Ld 

X 

U. 

300 °F TEMP. 


OPTIMIZATION VARIABLES 

INITIAL DESIGN 

FINAL DESIGN 

OBJECTIVE 

FUNCTION 

■ - 

NASTRAN BENDING 
MODULUS: G 1 1 
(ELEMENT 12/13) 

261.1580 

281.0925 

DESIGN 

VARIABLES 

PLY 

ANGLES 

(DEGREES) 

MATERIAL SYSTEM 0 

n 

0.0000 

0.0000 

2 

18.0000 

0.0040 

3 

36.0000 

- 0.1720 

D 

54.0000 

21.2350 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

D 

NOT USED 

NOT USED 

CONSTRAINT 

VARIABLES 

STRESS 

FAILURE 

CRITERIA 

F - 1.0 - MDElE 

NO FALURE : F < 1 
INCIPIENT : F « 1 
FALURE : F > 1 

MATERIAL SYSTEM § 

D 

0.0064 

0.0065 

2 

0.0243 

0.0031 

3 

0.0535 

0.0015 

B 

0.0788 

0.044 1 

5 

NOT USED 

NOT USED 

6 

NOT USED 

NOT USED 

B 

NOT USED 

NOT USED 
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